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ABSTRACT 

We revisit the distant future of the Sun and the solar system, based on stellar models 
computed with a thoroughly tested evolution code. For the solar giant stages, mass- 
loss by the cool (but not dust-driven) wind is considered in detail. Using the new and 
well-calibrated mass-loss formula of Schróder & Cuntz (2005, 2007), we find that the 
mass lost by the Sun as an RGB giant (0.332 Mo, 7.59 Gy from now) potentially gives 
planet Earth a significant orbital expansion, inversely proportional to the remaining 
solar mass. 

According to these solar evolution models, the closest encounter of planet Earth 
with the solar cool giant photosphere will occur during the tip-RGB phase. During 
this critical episode, for each time-step of the evolution model, we consider the loss 
of orbital angular momentum suffered by planet Earth from tidal interaction with 
the giant Sun, as well as dynamical drag in the lower chromosphere. As a result of 
this, we find that planet Earth will not be able to escape engulfment, despite the 
positive effect of solar mass-loss. In order to survive the solar tip-RGB phase, any 
hypothetical planet would require a present-day minimum orbital radius of about 1.15 
AU. The latter result may help to estimate the chances of finding planets around 
White Dwarfs. 

Furthermore, our solar evolution models with detailed mass-loss description pre- 
dict that the resulting tip-AGB giant will not reach its tip-RGB size. Compared to 
other solar evolution models, the main reason is the more significant amount of mass 
lost already in the RGB phase of the Sun. Hence, the tip-AGB luminosity will come 
short of driving a final, dust-driven superwind, and there will be no regular solar plan- 
etary nebula (PN). The tip-AGB is marked by a last thermal pulse and the final mass 
loss of the giant may produce a circumstellar (CS) shell similar to, but rather smaller 
than, that of the peculiar PN IC 2149 with an estimated total CS shell mass of just a 
few hundredths of a solar mass. 
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1 INTRODUCTION Apps 2001 (hereafter SSA)), and has been discussed very 


tly by Laughlin (2007). 
Climate change and global warming may have drastic effects pecently by Lamphun ) 


on the human race in the near future, over human time- 
scales of decades or centuries. However, it is also of interest, 
and of relevance to the far future of all living species, to 
consider the much longer-term effects of the gradual heating 
of the Earth by a more luminous Sun as it evolves towards 
its final stage as a white dwarf star. This topic has been 
explored on several occasions (e.g. Sackmann, Boothroyd & 
Kraemer 1993, Rybicki & Denis 2001, Schröder, Smith & 
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Theoretical models of solar evolution tell us that the 
Sun started on the zero-age main sequence (ZAMS) with 
a luminosity only about 70% of its current value, and it 
has been a long-standing puzzle that the Earth seems none 
the less to have maintained a roughly constant temperature 
over its life-time, in contrast to what an atmosphere-free 
model of irradiation would predict. Part of the explanation 
may be that the early atmosphere, rich in CO2 that was 
subsequently locked up in carbonates, kept the temperature 
up by a greenhouse effect which decreased in effectiveness 
at just the right rate to compensate for the increasing solar 
flux. The rôle of clouds, and their interaction with galactic 
cosmic rays (CR), may also be important: there is now some 
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evidence (Svensmark 2007; but see Harrison et al. 2007 and 
Priest et al. 2007) that cosmic rays encourage cloud cover 
at low altitudes, so that a higher CR flux would lead to a 
higher albedo and lower surface temperature. The stronger 
solar wind from the young Sun would have excluded galactic 
cosmic rays, so cloud cover on the early Earth may have been 
less than now, allowing the full effect of the solar flux to be 
felt. 

What of the future? Although the Earth’s atmosphere 
may not be able to respond adequately on a short time- 
scale to the increased greenhouse effect of carbon dioxide 
and methane released into the atmosphere by human activ- 
ity, there is still the possibility, represented by James Love- 
lock’s Gaia hypothesis (Lovelock 1979, 1988, 2006), that the 
biosphere may on a longer time-scale be able to adjust it- 
self to maintain life. Some doubt has been cast on that view 
by recent calculations (Scaife, private communication, 2007; 
for details, see e.g. Cox et al. 2004, Betts et al. 2004) which 
suggest that, on the century timescale, the inclusion of bio- 
spheric processes in climate models actually leads to an in- 
crease in carbon dioxide emissions, partly through a feed- 
back that starts to dominate as vegetation dies back. In any 
case, it is clear that the time will come when the increasing 
solar flux will raise the mean temperature of the Earth to a 
level that not even biological or other feedback mechanisms 
can prevent. There will certainly be a point at which life is 
no longer sustainable, and we shall discuss this further in 
Section B] 

After that, the fate of the Earth is of interest mainly 
insofar as it tells us what we might expect to see in systems 
that we observe now at a more advanced stage of evolution. 
We expect the Sun to end up as a white dwarf — do we expect 
there to be any planets around it, and in particular do we 
expect any small rocky planets like the Earth? 

The question of whether the Earth survives has proved 
somewhat tricky to determine, with some authors arguing 
that the Earth survives (e.g. SSA) and others (e.g. Sack- 
mann et al. 1993) claiming that even Venus survives, while 
general textbooks (e.g. Prialnik 2000, p.10) tend to say that 
the Earth is engulfed. A simple model (e.g. SSA), ignoring 
mass loss from the Sun, shows clearly that all the planets 
out to and including Mars are engulfed, either at the red 
giant branch (RGB) phase — Mercury and Venus - or at the 
later asymptotic giant branch (AGB) phase - the Earth and 
Mars. However, the Sun loses a significant amount of mass 
during its giant branch evolution, and that has the effect 
that the planetary orbits expand, and some of them keep 
ahead of the advancing solar photosphere. The effect is en- 
hanced by the fact (SSA) that when mass loss is included 
the solar radius at the tip of the AGB is comparable to that 
at the tip of the RGB, instead of being much larger; Mars 
certainly survives, and it appears (SSA) that the Earth does 
also. 

The crucial question here is: what is the rate of mass 
loss in real stars? Ultimately this must be determined from 
observations, but in practice these must be represented 
by some empirical formula. Most people use the classical 
Reimers’ formula (Reimers 1975, 1977), but there is consid- 
erable uncertainty in the value to be used for his parameter 
1), and different values are needed to reproduce the observa- 
tions in different parameter regimes. In our own calculations 
(SSA) we used a modification of the Reimers' formula, which 


has since been further improved and calibrated rather care- 
fully against observation, so that we believe that it is cur- 
rently the best available representation of mass loss from 
stars with non-dusty winds (Schróder & Cuntz 2005, 2007 
— see Section [2] where we explore the consequences of this 
improved mass-loss formulation). 

However, although we have considerably reduced the 
uncertainties in the mass-loss rate, there is another factor 
that works against the favourable effects of mass loss: tidal 
interactions. Expansion of the Sun will cause it to slow its 
rotation, and even simple conservation of angular momen- 
tum predicts that by the time the radius has reached some 
250 times its present value (cf. Table[I) the rotation period 
of the Sun will have increased to several thousand years in- 
Stead of its present value of under a month; effects of mag- 
netic braking will lengthen this period even more. This is so 
much longer than the orbital period of the Earth, even in 
its expanded orbit, that the tidal bulge raised on the Sun's 
surface by the Earth will pull the Earth back in its orbit, 
causing it to spiral inwards. 

This effect was considered by Rybicki & Denis (2001), 
who argued that Venus was probably engulfed, but that the 
Earth might survive. An earlier paper by Rasio et al. (1996) 
also considered tidal effects and concluded on the contrary 
that the Earth would probably be engulfed. However, the 
Rybicki & Denis calculations were based on combining an- 
alytic representations of evolution models (of Hurley, Pols 
& Tout 2000) with the original Reimers’ mass-loss formula 
rather than on full solar evolution calculations with a well- 
calibrated mass-loss formulation. The Rasio et al. paper also 
employed the original Reimers' formula, and both papers use 
somewhat different treatments of tidal drag. We have there- 
fore re-considered this problem in detail, with our own evolu- 
tionary calculations and an improved mass-loss description 
as the basis; full details are given in Sections [2] and [4] 


2 SOLAR EVOLUTION MODEL WITH MASS 
LOSS 


In order to describe the long-term solar evolution, we use the 
Eggleton evolution code (Eggleton 1971, 1972, 1973) in the 
version described by Pols et al. (1995, 1998), which has up- 
dated opacities and an improved equation of state. Among 
other desirable characteristics, his code uses a self-adapting 
mesh and a V-based prescription of “overshooting”, which 
has been well-tested and calibrated with giant stars in eclips- 
ing binaries (for details, see Schróder et al. 1997, Pols et al. 
1997, Schróder 1998). Because of the low mass and a non- 
convective core, solar evolution models are, however, not 
subject to any MS (main sequence) core-overshooting. In 
use, the code is very fast, and mass-loss is accepted simply 
as an outer boundary condition. 

As already pointed out by VandenBerg (1991), evolu- 
tion codes have the tendency to produce, with their most 
evolved models, effective temperatures that are slightly 
higher than the empirically determined values. The reason 
lies, probably, in an inadequacy of both low-temperature 
opacities and mixing-length theory at low gravity. With the 
latter, we should expect a reduced efficiency of the con- 
vective energy transport for very low gravity, because the 
largest eddies are cut out once the ratio of eddy-size to stellar 


radius has increased too much with g~*. Hence, as described 
by Schröder, Winters & Sedlmayr (1999), our mixing-length 
parameter, normally a = 2.0 for logg < 1.94, receives a 
small adjustment in the form of a gradual reduction for su- 
pergiant models, reaching a = 1.67 at logg = 0.0. With 
this economical adjustment, our evolution models now give 
a better match to empirically determined effective tempera- 
tures of very evolved late-type giants and supergiants, such 
as a’ Her (see Schröder & Cuntz 2007, Fig. 4 in particular), 
and even later stages of stellar evolution (Dyck et al. 1996, 
and van Belle et al. 1996, 1997). 

The evolution model of the Sun presented here uses an 
opacity grid that matches the empirical solar metallicity of 
Anders & Grevesse (1989), Z — 0.0188, derived from atmo- 
spheric models with simple 1D radiative transfer — an ap- 
proach consistent with our evolution models. Together with 
X = 0.700 and Y = 0.2812, there is a good match with 
present-day solar properties derived in the same way (see 
Pols et al. 1995). We note that the use of 3D-hydrodynamic 
modelling of stellar atmospheres and their radiative transfer 
may lead to a significantly lower solar abundance scale (e.g., 
Asplund, Grevesse & Sauval 2005, who quote Z — 0.0122), 
but these lower values are still being debated, and create 
some problems with helioseismology. Of course, using lower 
metallicities with an evolution code always results in more 
compact and hotter stellar models. Hence, if we used a lower 
Z our code would plainly fail to reproduce the present-day 
Sun, and the reliability of more evolved models with lower 
Z must therefore also be seriously doubted. 

The resulting solar evolution model suggests an age of 
the present-day MS Sun of 4.58 Gy (+0.05 Gy), counted from 
its zero-age MS start model, which is well within the range 
of commonly accepted values for the real age of the Sun and 
the solar system (e.g. Sackmann et al. 1993). Our model also 
confirms some well-established facts: (1) The MS-Sun has 
already undergone significant changes, i.e., the present solar 
luminosity L exceeds the zero-age value by 0.30 Lo, and the 
zero-age solar radius R was 11% smaller than the present 
value. (2) There was an increase of effective temperature To 
from, according to our model, 5596 K to 5774 K (+5 K). (3) 
The present Sun is increasing its average luminosity at a rate 
of 1% in every 110 million years, or 10% over the next billion 
years. All this is completely consistent with established solar 
models like the one of Gough (1981). 

Certainly, the solar MS-changes and their consequences 
for Earth are extremely slow, compared to the current cli- 
mate change driven by human factors. Nevertheless, solar 
evolution will force global warming upon Earth already in 
the “near” MS future of the Sun, long before the Sun starts 
its evolution as a giant star (see our discussion of the hab- 
itable zone in Section B). 

At an age of 7.13 Gy, the Sun will have reached its high- 
est Ter of 5820 K, at a luminosity of 1.26 Lo. From then on, 
the evolving MS Sun will gradually become cooler, but its 
luminosity will continue to increase. At an age of 10.0 Gy, 
the solar effective temperature will be back at Teg = 5751 K, 
while L = 1.84 Lo, and the solar radius then will be 37% 
larger than today. Around that age, the evolution of the 
Sun will speed up, since the solar core will change from cen- 
tral hydrogen-burning to hydrogen shell-burning and start 
to contract. In response, the outer layers will expand, and 
the Sun will start climbing up the RGB (the “red” or "first 
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Table 1. Main physical properties of characteristic solar models 


Phase Age/Gy  L/Lo9  To.g/K R/Ro Msun/Mo 
ZAMS 0.00 0.70 5596 0.89 1.000 
present 4.58 1.00 5774 1.00 1.000 
MS:hottest 7.13 1.26 5820 1.11 1.000 
MS:final 10.00 1.84 5751 1.37 1.000 
RGB:tip 12.17 2730. 2602 256. 0.668 
ZA-He 12.17 53.7 4667 11.2 0.668 
AGB:tip 12.30 2090. 3200 149. 0.546 
AGB:tip-TP 12.30 4170. 3467 179. 0.544 


(note: 1.00 AU = 215 Rọ) 


giant branch" in the HRD) - at first very gradually, but 
then accelerating. At an age of 12.167 Gy, the Sun will have 
reached the tip of the RGB, with a maximum luminosity of 
2730 Lo. 

In order to quantify the mass-loss rate of the evolved, 
cool solar giant at each time-step, we use the new and well- 
calibrated mass-loss formula for ordinary cool winds (i.e., 
not driven by dust) of Schréder & Cuntz (2005, 2007). This 
relation is, essentially, an improved Reimers’ law, physically 
motivated by a consideration of global chromospheric prop- 
erties and wind energy requirements: 


3.5 
M. \4000K 4300 gx 
with n = 8x107'*Mo y^ !, go = solar surface gravitational 
acceleration, and L., R+, and M. in solar units. 

This relation was initially calibrated by Schréder & 
Cuntz (2005) with the total mass loss on the RGB, using 
the blue-end (i.e., the least massive) horizontal-branch (HB) 
stars of globular clusters with different metallicities. This 
method avoids the interfering problem of temporal mass- 
loss variations found with individual giant stars and leaves 
an uncertainty of the new 7-value of only 15%, just under 
the individual spread of RGB mass-loss required to explain 
the width of HBs. 

Later, Schroder & Cuntz (2007) tested their improved 
mass-loss relation with six nearby galactic giants and su- 
pergiants, in comparison with four other, frequently quoted 
mass-loss relations. All but one of the tested giants are AGB 
stars, which have (very different) well-established physical 
properties and empirical mass-loss rates, all by cool winds 
not driven by radiation-pressure on dust. Despite the afore- 
mentioned problem with the inherent time-variability of this 
individual-star-approach, the new relation (equation (1)) 
was confirmed to give the best representation of the cool, 
but not “dust-driven” stellar mass-loss: it was the only one 
that agreed within the uncertainties (i.e., within a factor of 
1.5 to 2) with the empirical mass-loss rates of all giants. 
Hence, since the future Sun will not reach the critical lu- 
minosity required by a “dust-driven” wind (see Section B), 
we here apply equation (1) to describe its AGB mass-loss as 
well as its RGB mass-loss. 

The exact mass-loss suffered by the future giant Sun 
has, of course, a general impact on the radius of the solar 
giant, since the reduced gravity allows for an even larger 
(and cooler) supergiant. The luminosity, however, is hardly 
affected because it is mostly set by the conditions in the con- 
tracting core and the hydrogen-burning shell. In total, our 
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solar evolution model yields a loss of 0.332 Mo by the time 
the tip-RGB is reached (for 7 = 8 x 10! ^ Moy" +). This is 
a little more than the 0.275 Mc obtained by Sackmann et 
al. (1993), who used a mass-loss prescription based on the 
original, simple Reimers' relation. Furthermore, our evolu- 
tion model predicts that at the very tip of the RGB, the 
Sun should reach R = 256 Ro = 1.2 AU (see Fig. [1), with 
L = 2730Lco and Tew = 2602 K. More details are given in 
Table [1] 

By comparison, a prescription of the (average) RGB 
mass-loss rate with n = 7 x 10^ Moy" !, near the lower 
error limit of the mass-loss calibration with HB stars, yields 
a solar model at the very tip of the RGB with R = 249 R6, 
L = 27A2Lo, Ter = 2650 K, and a total mass lost on the 
RGB of 0.268Mo. With n = 9 x 10^" Moy" !, on the other 
hand, the Sun would reach the very tip of the RGB with 
R = 256Ro, L = 2714Lo, Ter = 2605 K, and will have lost 
a total of 0.388M. While these slightly different possible 
outcomes of solar tip-RGB evolution — within the uncer- 
tainty of the mass-loss prescription — require further discus- 
sion, which we give in Section [4.3] the differences are too 
small to be obvious on the scale of Fig. 

With the reduced solar mass and, consequently, lower 
gravitational attraction, all planetary orbits — that of the 
Earth included — are bound to expand. This is simply 
a consequence of the conservation of angular momentum 
Ag = Mg vg rg, while the orbital radius (i.e. rg) adjusts 
to a new balance between centrifugal force and the reduced 
gravitational force of the Sun, caused by the reduced so- 
lar mass Mgun(t). Substituting vg = 4/GMsun(t)/rg in Ag 
yields rg ox A% /Msun(t). For this conservative case, we find 
that rg is 1.50 AU for the case 7 = 8 x 10" May" !. For 
the smaller (7 x 10714) and larger (9 x 10714) values of n, 
we find, respectively, rg = 1.37 AU and rg = 1.63 AU, so 
in all cases the orbital radius is comfortably more than the 
solar radius, when angular momentum is conserved. 

Section provides a treatment of the more realistic 
case, in which angular momentum is not conserved. We have 
taken great care in determining the mass-loss and other pa- 
rameters for our models, because the best possible models of 
the evolution of solar mass and radius through the tip-RGB 
phase are required to provide reliable results. 

The significant solar RGB mass loss will also shape the 
later solar AGB evolution. Compared with models without 
mass loss, the AGB Sun will not become as large and lu- 
minous, and will be shorter-lived, because it lacks envelope 
mass for the core and its burning shells to “eat” into. In fact, 
the solar tip-AGB radius (149 Re) will never reach that of 
the tip-RGB (see Fig. [), and AGB thermal pulses are no 
threat to any planet which would have survived the tip- 
RGB. Our evolution code resolved only the two final and 
most dramatic thermal pulses (cf. Section [5]. 

The regular tip-AGB luminosity of 2090 Lo will not ex- 
ceed the tip-RGB value, either. Hence, as will be discussed 
in Section [5] the tip-AGB Sun will not develop a sustained 
dust-driven superwind but will stay short of the critical lu- 
minosity required by dust-driven winds (see Schróder et al. 
1999). The very tip of the AGB coincides with a thermal 
pulse (TP), after which the giant briefly reaches a peak lu- 
minosity of 4170 Lo, but at a higher Ter = 3467 K than on 
the RGB (see Table [1] and Section [5), keeping the radius 
down to 179 Ro. Again, the best possible treatment of all 
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age 
Figure 1. Solar radius evolution during the RGB and AGB 
phases. Included for comparison (dashed curve) is the potential 
orbital radius of planet Earth, taking account of solar mass loss 
but neglecting any loss of orbital angular momentum. The labels 
on the curve for the solar radius show the mass of the Sun in 
units of its present-day mass. 


prior mass loss from the giant Sun is essential for modelling 
this phase reliably. 


3 EVOLUTION OF THE HABITABLE ZONE 


The Earth currently sits in the ‘habitable zone’ in the so- 
lar system, that is, the region in which conditions on the 
Earth — in particular the average planetary temperature — 
are favourable for life. There are various precise definitions 
of ‘habitability’ in the literature, and a useful overview of 
habitable zones in the wider context of extrasolar planetary 
systems is given by Franck et al. (2002). For the current pa- 
per, a convenient definition is that a planet is habitable if 
the conditions on it allow the presence of liquid water on its 
surface. This may allow extremes of temperature that would 
make life uncomfortable if not impossible for humans, but 
the argument is that life of any kind (at least any kind we 
know about at present) requires water at some stage in its 
life cycle. We shall adopt that definition in this paper, but 
note that even with that apparently simple definition it is 
not straightforward to calculate the width of the habitable 
zone. 

It may be instructive to begin with a calculation of the 
mean planetary temperature in terms of a spherical black 
body, by assuming that the planetary body absorbs the so- 
lar flux intercepted by its (circular) cross-sectional area and 
re-emits it spherically symmetrically at a black body tem- 
perature T'. Then (cf. SSA) T' is given by 


1/2 
T = a-a (25) Tort 
1/2 1/2 
R 1AU 
E _ Ay A 
— 0.0682 (1— A) (£) (=) Ter — (2) 


where D is the distance of the body from the centre of the 
Sun, R is the radius of the Sun, A is the Bond albedo of the 
Earth and Teg is the effective temperature of the Sun. On 
that basis, taking Tog = 5774 K and R = Ro (Table[1], and 
A = 0.3 (Kandel & Viollier 2005), we find T(1 AU) = 255 K. 


But the actual mean temperature of the Earth at present is 
33 K warmer, at T = 288 K. This demonstrates the warming 
effect of our atmosphere, which becomes significantly more 
important with higher temperature (see below). 

In fact, there are various complex, partly antagonistic 
atmospheric feedback mechanisms (for example, the green- 
house effect, the variation of planetary albedo with the pres- 
ence of clouds, snow and ice, and the carbonate-silicate cycle 
which determines the amount of carbon dioxide in the at- 
mosphere) that act to change the surface temperature from 
what it would be in the absence of an atmosphere. These 
mechanisms have been carefully discussed by Kasting, Whit- 
mire & Reynolds (1993), who conclude that a conservative 
estimate of the current habitable zone (HZ) stretches from 
0.95 AU to 1.37 AU. We shall adopt their result for the lim- 
ited purposes of this paper. It can be adjusted in a simple- 
minded way to allow for the evolution of the Sun by scaling 
the inner and outer HZ radii ryz,i, rnz,o with the changing 
solar luminosity Lsgun(t): raz x 4/ Lsux(t). In this way, the 
respective critical values of solar irradiance derived by Kast- 
ing et al. (1993) for the inner and outer edge of the HZ are 
maintained. 

Certainly, with the 10% increase of solar luminosity over 
the next 1Gy (see previous section), it is clear that Earth 
will come to leave the HZ already in about a billion years 
time, since the inner (hot side) boundary will then cross 
1AU. By the time the Sun comes to leave the main se- 
quence, around an age of 10 Gy (Table[I), our simple model 
predicts that the HZ will have moved out to the range 1.29 
to 1.86 AU. The Sun will have lost very little mass by that 
time, so the Earth's orbital radius will still be about 1 AU — 
left far behind by the HZ, which will instead be enveloping 
the orbit of Mars. 

By the time the Sun reaches the tip of the RGB, at 
12.17 Gy, the Earth’s orbital radius will only have expanded 
to at most 1.5 AU, but the habitable zone will have a range 
of 49.4 to 71.4 AU, reaching well into the Kuiper Belt! The 
positions of the HZ boundaries are not as well determined 
as these numbers suggest, because in reality the scaling for 
the boundaries of the HZ almost certainly depends also on 
how clouds are affected by changes in the solar irradiance. 
These effects are complex and uncertain (cf. Kasting 1988), 
and may increase or decrease the speed at which the HZ 
drifts outwards. But none the less it seems clear that the 
HZ will move out past the Earth long before the Sun has 
expanded very much, even if the figure of one billion years 
is a rather rough estimate of how long we have before the 
Earth is uninhabitable. 

In other planetary systems around solar-type stars, con- 
ditions may be different, and it may even be possible for life 
to start during a star’s post-main-sequence evolution, if a 
planet exists at a suitable distance from the star. This pos- 
sibility is discussed by Lopez, Schneider & Danchi (2005), 
who also give a general discussion of the evolution of habit- 
able zones with time. However, they use the evolution mod- 
els of Maeder & Meynet (1988), which do not agree as well 
as ours with the colours and observed Tem’s of the red gi- 
ants in star clusters (see, e.g., illustrations given by Meynet 
et al. 1993), and which predict a very different behaviour for 
the solar radius; so their results are not directly comparable 
with ours. 

What will happen on the Earth itself? Ignoring for the 
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moment the short-timescale (decades to centuries) problems 
currently being introduced by climate change, we may ex- 
pect to have about one billion years of time before the solar 
flux has increased by the critical 1096 mentioned earlier. At 
that point, neglecting the effects of solar irradiance changes 
on the cloud cover, the water vapour content of the atmo- 
sphere will increase substantially and the oceans will start 
to evaporate (Kasting 1988). An initially moist greenhouse 
effect (Laughlin 2007) will cause runaway evaporation until 
the oceans have boiled dry. With so much water vapour in 
the atmosphere, some of it will make its way into the strato- 
sphere. There, solar UV will dissociate the water molecules 
into OH and free atomic hydrogen, which will gradually es- 
cape, until most of the atmospheric water vapour has been 
lost. The subsequent dry greenhouse phase will raise the sur- 
face temperature significantly faster than would be expected 
from our very simple black-body assumption, and the ulti- 
mate fate of the Earth, if it survived at all as a separate 
body (cf. Section 4), would be to become a molten remnant. 


4 THE INNER PLANETARY SYSTEM 
DURING TIP-RGB EVOLUTION 


After 12 Gy of slow solar evolution, the final ascent of the 
RGB will be relatively fast. The solar radius will sweep 
through the inner planetary system within only 5 million 
years, by which time the evolved solar giant will have 
reached the tip of the RGB and then entered its brief (130 
million year) He-burning phase. The giant will first come to 
exceed the orbital size of Mercury, then Venus. By the time 
it approaches Earth, the solar mass-loss rate will reach up to 
2.5 x 107" Mc y^! and lead to some orbital expansion (see 
Section[2). But the extreme proximity of the orbiting planet 
to the solar photosphere requires the consideration of two 
effects, which both lead to angular momentum loss and a 
fatal decrease of the orbital radius of planet Earth. 


4.1 Tidal interaction 


For the highly evolved giant Sun, we may safely assume 
(cf. Section [I] that it has essentially ceased to rotate, after 
nearly 2 billion years of post-MS magnetic braking acting on 
the hugely expanded, cool RGB giant. Consequently, any 
tidal interaction with an orbiting object will result in its 
suffering a continuous drag by the slightly retarded tidal 
bulges of the giant solar photosphere. 

As shown in Section[2] the orbital radius of planet Earth 
rp depends on the angular momentum squared, by the equa- 
tion 


Ax (t) 


~ MZ G Msat) (3) 


TE 
Hence, the terrestrial orbit reacts quite sensitively to any 
loss of angular momentum, by shrinking. 

The retardation of the tidal bulges of the solar photo- 
sphere will be caused by tidal friction in the outer convec- 
tive envelope of the RGB Sun. This physical process was 
analyzed, solved and applied by J.-P. Zahn (1977, 1989, and 
other work referred to therein), and successfully tested with 
the synchronization and circularization of binary star orbits 
by Verbunt & Phinney (1995). In a convective envelope, the 
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main contribution to tidal friction comes from the retarda- 
tion of the equilibrium tide by interaction with convective 
motions. For a circular orbit, the resulting torque I’ exerted 
on planet Earth by the retarded solar tidal bulges is given 
by (Zahn 1977; Zahn 1989, Eq.11): 


Ap Rsun 
tr TE 


6 

r=6 Mss Ran ( ) (9 — w). (4) 
Here, the angular velocity of the solar rotation is sup- 
posed to be Q = 0, while that of the orbiting Earth, 
w(t) = 2n/ Pz(t) = A ?(t) ME&(GMsus(t))?, will vary both 
with the decreasing angular momentum A(t) (— 2.67 x 
10€ kgm?s^! at present) and with the solar mass in the 
final solar RGB stages. The exerted torque scales with 
the square of the (slowly increasing) mass ratio q(t) — 
Mr/Msun(t) (= 3.005 x 10 9 at present), because q 
determines the magnitude of the tidal bulges. t;(t) — 
(Msux (£) Ran (t) /Lsun(t))'/? ~ O(1y) is the convective fric- 
tion time (Zahn 1989, Eq.7), and the coefficient A2 depends 
on the properties of the convective envelope. For a fully con- 
vective envelope (Zahn 1989, Eq.15), with a tidal period ~% 
O(1y), comparable to 2£;, we may use Ao z 0.019 atl? a 
0.038 (with a convection parameter of our tip-RGB solar 
model of a ~ 1.7). This coefficient appears to be the main 
source of uncertainty (see Section [4.3), because it is related 
to the simplifications of the mixing length theory (MLT). 

With the properties of the tip-RGB Sun, a typical value 
of the tidal drag acting on planet Earth is T = dA/dt = 
—3.3 x 1079 kg m? s^?, which gives a typical orbital angular 
momentum decay time of r =| A/T |= 2.6 x 10° y. This is 
comparable to the time spent by the Sun near the tip-RGB; 
since a loss of only ~ 10% of the angular momentum will 
be sufficient to reduce the orbital radius (by 20%) to lower 
it into the solar giant photosphere, this order-of-magnitude 
calculation illustrates clearly that tidal interaction is crucial. 
Its full consideration requires a timestep-by-timestep com- 
putation of the loss of orbital angular momentum; at each 
time-step of the solar evolution calculation, we use equa- 
tion (4), together with the radii and masses of our solar 
evolution model, to compute the change in angular momen- 
tum, and then use equation to compute the change in 
the orbital radius, and hence the new orbital period of the 
Earth. Section presents the result, which also takes into 
account the relatively small additional angular momentum 
losses by dynamical drag, as discussed in the next section. 


4.2 Dynamical friction in the lower chromosphere 


A further source of angular momentum loss by drag is dy- 
namical friction, from which any object suffers in a fairly 
close orbit, by its supersonic motion through the gas of 
the then very extended, cool solar giant chromosphere. In a 
different context, dynamical drag exerted by a giant atmo- 
sphere has already been considered by Livio & Soker (1984). 
But the specific problem here is to find an adequate de- 
scription of the density structure of the future cool solar gi- 
ant. Fortunately, as it turns out (see below), dynamical drag 
will play only a minor róle, very near the solar giant photo- 
sphere, and the total angular momentum loss is dominated 
by the tidal interaction described above. An approximate 
treatment of the drag is therefore adequate, and we use the 
recent study by Ostriker (1999). 


In the case of supersonic motion (with a Mach num- 
be] of the order of 2 to 3) in a gaseous medium, dynamical 
friction consists in about equal shares of the collisionless, 
gravitational interaction with its wake and of the friction it- 
self. In her study, Ostriker (1999, Fig. 3) finds that the drag 
force exerted on the object in motion is 


Fy = \a4mp(GMp/cs)” (5) 


where A, is of the order of 1 to 3. The numerical simulations 
made by Sánchez-Salcedo & Brandenburg (2001) are in gen- 
eral agreement with the results of Ostriker (1999). Here c; 
is the speed of sound, which in a stellar chromosphere is 
about 8kms !, and p is the gas density (SI units). The lat- 
ter quantity is the largest source of uncertainty, as we can 
only make guesses (see below) as to what the gas density 
in the lower giant solar chromosphere will be. The angular 
momentum loss resulting from this drag is simply 


dA/dt = —Farg, (6) 


and the corresponding life-time of the orbital angular mo- 
mentum is T = A/ | dA/dt |, as above. 

For the lower chromosphere of the K supergiant ¢ Aur, 
employing an analysis of the additional line absorption in 
the spectrum of a hot companion in chromospheric eclipse, 
Schróder, Griffin & Griffin (1990) found an average hydrogen 
particle density of 7 x 10™ cm? at a height of 2 x 108 km. 
Alternatively, we may simply assume that the density of the 
lower solar chromosphere scales with gravity g, which will be 
lower by 4.7 orders of magnitude on the tip-RGB, while the 
density scale-height scales with g^! (as observations of cool 
giant chromospheres seem to indicate, see Schróder 1990). 
The chromospheric models of both Lemaire et al. (1981) and 
Maltby et al. (1986) suggest particle densities of the order of 
10!" cem"? at a height of 100 km, and a scale height of that 
order for the present, low solar chromosphere. Scaled to tip- 
RGB gravity, that would correspond to a particle density 
of 2 x 10? cm^?, or p ~ 4 x 107? kgm ?, at a height of 
5 x 106 km (0.03 AU), and a density scale height of that 
same value. 

For the computation of the orbital angular momentum 
loss of the Earth, presented below (see Figures 2]and[3), we 
apply the latter, rather higher values of the future chromo- 
spheric gas density, together with the (also more pessimistic) 
assumption of Ag = 3 (using c; = 8km s^ 1). The typical an- 
gular momentum decay-time by dynamical friction in the 
low (h © 0.03 AU) chromosphere of the tip-RGB solar giant 
is 14 million years — significantly longer than that for tidal 
interaction. Hence, this illustrates that dynamical friction is 
of interest only in the lowest chromospheric layers, adding 
there just a little to the drag exerted by tidal interaction. 
None the less, we include it, using equations (5) and (6) to 
calculate the additional angular momentum change to be 
included in equation (8). 


1 Note that vg c Msun(t), and so the Mach number is somewhat 
lower than would be expected from the present orbital velocity of 
the Earth of about 30 km s-!. 
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Figure 2. The final 4 million years of solar evolution before the 

tip-RGB, showing the radii of the Sun and of the orbit of planet 

Earth (dashed curve) — taking account of angular momentum 

losses by tidal interaction and supersonic drag. The labels on the 

solar radius track give values of Mgun(t)/Mo, as in Figure[1] 


4.3 “Doomsday” confirmed 


As explained in the previous two sections, we use equations 
to (6) to compute, at each time-step of our evolutionary 
calculation, a detailed description of the orbital evolution for 
planet Earth in the critical tip- RGB phase of the Sun under 
the influence of tidal interaction and dynamical drag. The 
resulting evolution both of the orbital radius of the Earth 
and of the radius of the solar giant is shown in Fig.|2] This 
shows that, despite the reduced gravity from a less massive 
tip-RGB Sun, the orbit of the Earth will hardly ever come to 
exceed 1 AU by a significant amount. The potential orbital 
growth given by the reduced solar mass is mostly balanced 
and, eventually, overcome by the effects of tidal interaction. 
Near the very end, supersonic drag also becomes a significant 
source of angular momentum loss. 

As shown by Fig.[2] engulfment and loss of planet Earth 
will take place just before the Sun reaches the tip of the 
RGB, 7.59 Gy (+0.05 Gy) from now. According to our calcu- 
lation, it occurs when the RGB Sun has still another 0.25 AU 
to grow, about 500,000 years before the tip-RGB. Of course, 
Mercury and Venus will already have suffered the same fate 
as Earth some time before — respectively, 3.8 and 1.0 million 
years earlier. 

As mentioned in the introduction, a similar calculation 
was already carried out in the context of extra-solar plan- 
ets by Rasio et al. (1996), who basically came to the same 
conclusions; their fig. 2 is reminiscent of ours. They also em- 
ployed the orbital decay rate predicted by Zahn's theory, but 
their solar evolution model used the old Reimers mass-loss 
relation, and they did not make any adjustments to match 
the effective temperatures found empirically at the tip of the 
giant branches (see Section [2] . 

Do the remaining uncertainties allow the possibility for 
Earth to escape the “doomsday” scenario? As far as the 
mass-loss alone is concerned, this seems unlikely: according 
to the study of HB stars in globular clusters by Schróder 
& Cuntz (2005), is remarkably well constrained and can- 
not exceed 9 x 10^ Ma y^ !, or the total RGB mass-loss 
would become so large that the tip-RGB star would miss He- 
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ignition and not reach the horizontal branch at all. And the 
full width of the HB towards lower Tug is achieved already 
with an 5 of 7 x 10^ Mo y~". Furthermore, the benefit of 
larger orbits with a reduced solar mass is to some extent 
compensated for by a larger solar giant. 

Dynamical drag does not become important until the 
planet is already very near the photosphere, i.e., after tidal 
drag has already lowered the orbit. Hence, the most signif- 
icant uncertainty here comes from the scaling of the tidal 
friction coefficient A» (of Zahn, 1989). For this reason, we 
computed several alternative cases, and from these we find: 

(1) With the mass-loss rate unchanged, the value of A2 
would have to be significantly smaller for an escape from 
the “doomsday” scenario, i.e., less than 0.013, instead of 
our adopted value of 0.038. But Zahn's scaling of A» has 
been empirically confirmed within a factor of 2, if not bet- 
ter (see Verbunt & Phinney, 1995). Very recently, realistic 
3D simulations of the solar convection have also resulted in 
an effective viscosity which matches that of Zahn's prescrip- 
tion surprisingly well (Penev et al. 2007). And Rybicki & 
Denis (2001), by comparison, used a value (2 = 0.05 in 
the notation of their very similar calculation of tidal angu- 
lar momentum loss) which is entirely consistent with Zahn's 
scaling of A» for a convection parameter of a = 2. 

(2) We then considered solar evolution models with a 
reasonably larger mass-loss rate (7 = 9 x 10 Mo y^!) in 
combination with tidal friction coefficients of 1/1, 2/3 and 
1/2 of the one given by Zahn. In each of these cases, planet 
Earth would not be able to escape doomsday but would 
face a delayed engulfment by the supergiant Sun — 470,000, 
230,000 and 80,000 years before the tip-RGB is reached, 
respectively. 

(3) Finally, we checked the outcome for a reasonably 
lower mass-loss rate (1 = 7 x 10! Mo y^!) in combination 
with the same tidal friction coefficients as above. The en- 
gulfment would then happen rather earlier than with more 
mass-loss — 540,000, 380,000 and 270,000 years before the 
tip-AGB is reached. 

These computations confirm that reducing the solar 
mass enlarges the planetary orbit more than the tip-RGB 
solar radius, so that the best way to avoid the doomsday 
scenario would be to have as high a mass-loss rate as possi- 
ble. However, we believe that the value of ņ in case (2) above 
already is as high as it can be without violating agreement 
of evolved models with observations, and that the smallest 
value used there for the tidal friction coefficient is also at 
the limits of what is allowed by observational constraints. 
The only possible escape would be if our solar giant models 
were too cool (by over 100 K in case 2), and therefore larger 
than the real Sun will be. Hence, to avoid engulfment by the 
tip-RGB Sun would require that all three parameters (7, A2 
and Tep) were at one edge of their uncertainty range, which 
seems improbable. Rather, our computations confirm, with 
reasonable certainty, the classical “doomsday” scenario. 


4.4 “Doomsday” avoidable? 


Even though this is an academic question, given the hos- 
tile conditions on the surface of a planet just missing this 
“doomsday” scenario, we may ask: what is the minimum ini- 
tial orbital radius of a planet in order for it to “survive”? 
Fig.[3]shows, by the same computation as carried out for Fig. 
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Figure 3. As Fig.[2] but for a planet with a present orbital radius 

of 1.15 AU. 


that an initial orbital radius of 1.15 AU is sufficient for 
any planet to pass the tip-RGB of a star with M; = 1.0 Mo. 
Since, as shown in Section[5] the tip-AGB Sun will not reach 
any similarly large extent again, such a planet will eventu- 
ally be orbiting a White Dwarf. 

A more general discussion of planetary survival during 
post-main-sequence evolution has been given by Villaver & 
Livio (2007), who suggest that an initial distance of at least 
3 AU is needed for the survival of a terrestrial-size planet 
when one also takes into account the possible evaporation 
of the planet by stellar heating. However, they use stellar 
models and mass-loss rates that have the maximum radius 
and mass loss occurring on the AGB. That has been the 
expected result for many years, but is quite different from 
what we find (Section [5] and Table 1) with the improved 
mass-loss formulation of Schróder & Cuntz (2005, 2007). 
Hence, Villaver & Livio's results may be unduly pessimistic. 

In any case, it is clear that terrestrial planets can survive 
if sufficiently far from their parent star. If it were possible 
to increase the orbital radius from its initial value, then an 
increase of only 896 of angular momentum should yield the 
pre-RGB orbital size required by planet Earth to escape 
engulfment. Is that conceivable? 

An ingenious scheme for doing so which, in the first 
place, could increase the time-scale for habitation by in- 
telligent life for the whole of the Sun's MS life-time, was 
proposed by Korycansky, Laughlin & Adams (2001). They 
pointed out that a suitable encounter of the Earth every 
6000 years or so with a body of large asteroidal mass could 
be arranged to move the orbit of the Earth outwards; Kuiper 
Belt objects might be the most suitable. The energy require- 
ments could be reduced by incorporating additional encoun- 
ters with Jupiter and/or Saturn. Although still very large 
by today's standards, the energy requirements remain small 
compared to those for interstellar travel. 

On the face of it, this scheme seems far-fetched, but 
Korycansky et al. (2001) show that it is in principle pos- 
sible, both technically and energetically, although currently 
somewhat beyond our technical capabilities; however, there 
is no immediate hurry to implement the scheme, which could 
await the development of the relevant technology. It would 
have the advantage of improving conditions for the whole 
biosphere, whereas any scheme for interplanetary ‘life rafts’ 
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Figure 4. Solar mass loss during the final 1 million years on the 
AGB will remain mainly of the order of 2 x 1077 Mo y- ! and not 
provide sufficient CS shell mass to form a regular PN. Only the 
last two TP’s (tip-AGB and post- AGB, see text) are resolved. 


that could move slowly outwards to maintain habitable con- 
ditions would, on cost and energy grounds, necessarily be 
confined to a small fraction of the human population — with 
all the political problems that that would produce — plus 
perhaps a tiny proportion of other species. None the less, 
the asteroidal fly-by scheme has its own problems, not least 
the danger of a benign close approach turning into a catas- 
trophic accidental collision, and possibly also triggering or- 
bital instability — cf. also Debes & Sigurdsson (2002). 


5 TIP-AGB SOLAR EVOLUTION 


The loss of 1/3 of the solar mass during the rise to the tip of 
the RGB will make a significant impact on the further evolu- 
tion as an AGB star. There is very little shell mass left, into 
which the two burning shells (H, followed by He) can ad- 
vance (on a radial mass scale). Hence, the C/O core cannot 
grow as much as with a conservative model without mass 
loss, and the whole core region will not contract as much, 
either. Consequently, the AGB luminosity, determined by 
the density and temperature in the H-burning shell, will not 
reach as high levels as in a conservative AGB model, and nei- 
ther will the AGB radius of the late future Sun (see Table 
fi). 

According to our evolution model, the regular tip-AGB 
evolution will be ended with a luminosity of only 2090 Lo, 
Tef = 3200 K, and R = 149 Re. The AGB mass-loss rate, 
according to the relation of Schróder & Cuntz (2005), will 
reach only 2.0 x 107’ Mc y^! (see Fig. Ø), since the lumi- 
nosity will not be sufficient to drive a dust-driven wind (see 
Schróder et al. 1999). Also, even if it did: only a little shell 
mass will have been left to lose after the RGB phase, only 
0.116Mo. 

Hence, for this non-dust-driven AGB solar mass-loss, 
we have adopted the same mass-loss description as above 
(equation (1)). This mass-loss, in combination with our solar 
evolution model, yields the following prediction: during the 
final 30,000 y on the very tip-AGB, which are crucial for any 
build-up of sufficient CS (circumstellar) material to form 
a PN, the solar giant will lose only 0.006 Mo. A further 


0.0015 Mo will be lost in just 1300 years right after a final 
thermal pulse (TP) on the tip-AGB. That marks the very 
end of AGB evolution, and it allows the solar supergiant 
briefly to reach a luminosity of 4170 Lọ and R = 179 Ro, 
with a mass-loss rate of 107° Mo y~', but with Tyg already 
increased to 3467 K. Again, there will be no involvement 
of a dust-driven wind. Since common PNe and their dusty 
CS envelopes reveal a dust-driven mass-loss history of more 
like 107° to 10^ Mc y^! during the final 30,000 years of 
tip-AGB evolution, we must conclude that the Sun will not 
form such a PN. 

Since a circumstellar shell of nearly 0.01 Mo will, nev- 
ertheless, be produced by the tip- AGB solar giant, a rather 
peculiar PN may be created by the emerging hot stellar core 
— it might be similar to IC 2149. Although most of the pe- 
culiar, strongly bi-polar PNe appear to stem from massive 
stars, this particular object has only a slim total mass of 
0.01 to 0.03 Mo, lacking a massive envelope — see Vázquez 
et al. 2002. Hence, these authors argue that this PN appears 
to be the product of a low-mass star with M; close to 1 Mo. 

A final mass of 0.0036 Mo is lost by the post-AGB 
star, which on its way to become a hot subdwarf undergoes 
at least one more TP. For the resulting solar WD (white 
dwarf), our evolution model yields a final mass of 0.5405 Mo. 


6 CONCLUSIONS 


We have applied an improved and well-tested mass-loss re- 
lation to RGB and AGB solar evolution models, using a 
well-tested evolution code. While the habitable zone in the 
inner solar system will already move outwards considerably 
in the next 5 billion years of solar MS evolution, marking 
the end of life on Earth, the most critical and fatal phase for 
the inner planetary system is bound to come with the final 
ascent of the Sun to the tip of the RGB. 

Considering in detail the loss of angular momentum by 
tidal interaction and dynamical drag in the lower chromo- 
sphere of the solar giant, we have been able to compare the 
evolution of the RGB solar radius with that of the orbit of 
planet Earth. Our computations reveal that planet Earth 
will be engulfed by the tip-RGB Sun, just half a million 
years before the Sun will have reached its largest radius of 
1.2 AU, and 1.0 (3.8) million years after Venus (and Mer- 
cury) have suffered the same fate. While solar mass loss 
alone would allow the orbital radius of planet Earth to grow 
sufficiently to avoid this “doomsday” scenario, it is mainly 
tidal interaction of the giant convective envelope with the 
closely orbiting planet which will lead to a fatal decrease of 
its orbital size. 

The loss of about 1/3 of the solar mass already on the 
RGB has significant consequences for the solar AGB evo- 
lution. The tip-AGB Sun will not qualify for an intense, 
dust-driven wind and, hence, will not produce a regular PN. 
Instead, an insubstantial circumstellar shell of just under 
1/100 Mc will result, and perhaps a peculiar PN similar to 
IC 2149. 
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